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Abstract 

The calculation presented in "A neoclassical calculation of rotation profiles and comparison with 
DIII-D measurements" by Stacey, Johnson, and Mandrekas, [Physics of Plasmas, 13, (2006)], 
contains several errors, including the neglect of the toroidal electric field, an unphysical expression 
for the electrostatic potential, and an unevaluated relation among its parameters. An alternative 
formulation is discussed. 
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The calculation presented in "A neoclassical calculation of rotation profiles and compar- 
ison with Dlll-D measurements" by Stacey, Johnson, and Mandrekas [l], contains several 
errors, including the neglect of the toroidal electric field, an unphysical expression for the 
electrostatic potential, and an unevaluated relation among its parameters. These inaccu- 
rate and neglected expressions have been part of the Georgia Tech Fusion Research Center 
(GT-FRC) model for over fifteen years 0, l^l and hamper its ability to interpret the results 
of experimental investigation. 

Reference [l] writes the toroidal equation of motion for species j as 
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While the toroidal electric field is retained symbolically in Equation (27) thereof, in prac- 
tice the term is unevaluated, as no Ohm's law (generalized or otherwise) appears in the 
presentation of the theory. Thus, its determination of the gyroviscous angular momentum 
transport frequency via the solution of Equation (31) thereof is cast in doubt, as the re- 
sulting force balance equation neglects a term omnipresent for inductive current drive. A 
consistent treatment of the toroidal electric field in a plasma of finite conductivity necessarily 
involves retention of the full electron equation of motion with finite electron mass. Following 

n 

Sturrock [4[], the conductivity tensor in the concentric circular flux surface approximation, 
~B/B = (0, b e , b^), may be written using coordinate rotation tensors as 
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given component conductivities an ~ n e e 2 /m e i/ ee , a±_ 
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( ?±(u e /b , ee ), for (unsigned) electron gyrofrequency u e and collision frequency u ee . Thus, 
the toroidal electric field induces currents in the poloidal plane as well as the toroidal direc- 



tion, 



j&</> j-C'tp -t^4> 



(bean, b^,be(a\\ — a±), b\a\\ + b 2 e a±) E$. With experimental 
determination of = +a\\V r Be from the equilibrium equation Vp = J x B and the sin- 
gle fluid velocity V = J2j n j m j^j/J2j n j m j^ one ma y determine E^ = / (b^a\\ + bga±). 
Note that the radial current (including terms neglected by Sturrock), J r = r ■ a ■ E + 
ay [dp e /n e edr + (VeB^ — V^Bg)], must be zero at equilibrium. 

Reference |l| writes the poloidal equation of motion for species j as 
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and takes the poloidal electrostatic field on a flux surface at r as given by 

ld$(r) 1 d n $°(r) 

E e = ^ = — —<S>°(r)(l + § c cos9 + $ s sm9) = ^ ($ s costf - $ c sin 9) , 

r o9 r 09 r 

(4) 

where $ is the electrostatic potential, indicating an expansion around $°(r) = — J drE® ^ 
for a last closed flux surface at r = a. The resulting evaluation of the flux surface unity, 
cosine, and sine moments of the electron poloidal equation of motion T e dn e /rd6 = —en e Eo 
(other terms are assumed negligible at equilibrium), defined by the expressions (A)u,c,s = 
§ d9{l, cos 9, sin 9}(1 + e cos 9)A/2tt, yields three equations which have only trivial solution. 
(The evaluation of these moments is best accomplished without recourse to the logarithmic 
derivative.) Specifically, for n e = n° e (l + n c e cos 9 + n s e sin 9), we have the equations 

U : e%n s e = e$°(e$ s + <$ s - <$ c ) , (5) 

C: %n s e = e$°(4$ s + 3en c e $ s - en s e ^ c )/A , (6) 

S: %n c e = e$°(4$ c + e<$ c -e<<l> s )/4 , (7) 

valid V e,nl,n s e . Independent variation thereof for finite (fixed) $° yields $ c,s = for this 
overdetermined system, thus the poloidal electrostatic field in the GT-FRC model, which 
fails to consider the unity moment equation and the 0(e) terms, must vanish. Failing to 
include the 0(e) terms indicates expressions applicable only on the magnetic axis r = 0, 
where e = t/Rq = for Rq ^ oo, yet the GT-FRC model is commonly used to address 
the physics near the edge of the confinement region 5j. These equations may be linearized 
and have the formal solution (1/$°, $ c , $ s ) = (0, 0, 0), which we interpret to mean exactly 
what it says, that they are solved when $°(r) = oo, displaying its unphysical definition. As 
non- vanishing $ c ' s are an integral part of the development of the GT-FRC model, appearing 
in both the poloidal and toroidal equations of motion, the validity of its conclusions is in 
jeopardy. 

Note in passing that from n e = ^2 ions Zini and dp e /d9 = — J2ions I '> which for two 
ion species Zi = 1 and z z ^ 1 give us the pair of equations 

n? = nf s + z z nT = -(pVpIW ~ {pVp>T , (8) 

the poloidal variations of the electron density, nl ,s , must vanish for a "pure" plasma consist- 
ing of a single (or single "effective") ion species but are supported in an "impure" plasma 
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consisting of multiple distinct ion species, when the electron and ion temperatures assume 
no poloidal variation as in the GT-FRC model. A single ion species plasma may support 
poloidal density asymmetries given by the poloidal dependence of the ion temperature, 
dlogrii/dO = —d\og(T e +Ti)/d9, which may also impact the multi-species treatment. The 
rightmost equality of Equations (|S]), providing an additional relation between n^ s and n c z ,s , 
is neglected in the GT-FRC model, which treats the ion and impurity density asymmetries 
as independent free parameters. 

An alternative treatment of the problem begins with the consideration of the poloidal 
line integral of the electric field, § dl ■ E = § rdOEg = at equilibrium. (Note that a 
non-vanishing radial electrostatic field without poloidal variation demands the existence of 
no poloidal electrostatic field, else the poloidal variation to the potential ruins the poloidal 



symmetry of the radial fie 
poloidal electrostatic field 



d.) Examining the expression of the leading contender for the 
6| evaluated from the equation of motion, 
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one may put it in the form Eg = E^[2ecos6 — (s 2 /2) cos 26], when the Shafranov shift is 
neglected, as in the concentric circular flux surface approximation of the GT-FRC model 
(and others), by requiring § dOEg = 0. Inserting that expression into the electron poloidal 
equation of motion and taking the flux surface Fourier moments yields three equations which 
have a nontrivial solution only when expanded to order 0(e 3 ), given by 
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thus the presence of a poloidal electrostatic field of that form should be accompanied by 
a potentially measurable shift in the electron density profile. Note that the derivation 
immediately preceding is slightly inconsistent, as the associated electrostatic potential $ oc 
2esin# — (e 2 /A) sin 26 is of the correct harmonic form — VL^ = —(d 2 /dZ 2 + d 2 /dR 2 )& = 



(rd/dr)/rdr + d 2 /r 2 d6 2 }§ = for axial geometry 



10| as implied by the neglect 



of the Shafranov shift, yet the flux surface average is done in toroida 



by taking the potential over to tokamak geometry is forthcoming 11] 



geometry. Resolution 
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